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1. Introduction 

The Universe is flat, homogeneous and isotropic, and has the primordial density perturba- 
tions with almost scale-invariant spectrum. 1 All these above may be well explained by a 
hypothetical inflationary stage of the Early Universe before it became hot. This stage can 
be realized within the Standard Model (SM) of particle physics extended by non-minimal 
coupling of the SM Higgs field H to gravity [1] 

^gravity = ~£-R + , (1.1) 

where M P = 1/^8itG n = 2.44 x 10 18 GeV is the reduced Planck mass. In [2] (see also 
[3]) we analyzed the history of the Universe within the Higgs inflation scenario in two 
theories, assuming their validity up to the Planck scale. The first one was the Standard 
Model and the second one was the z/MSM (an extension of the SM by 3 singlet fermions 

1 lt is well known for a long time, e.g.: « Hyduxufi o6jio, 03opno, ozpoMno c mpu3eeH0u u JIaeu,» B. Tpe- 
flnaKOBCKHH, TejieMaxuda (1766). 



- 1 - 



with masses below the electroweak scale, for a review see ref. [4]). One of our main goals 
in [2] was to determine the initial conditions for the hot Big Bang. In other words, we 
addressed the question whether the quantum numbers, effectively conserved in the SM 
(such as B — L, where B and L are baryon and lepton numbers respectively) and in the 
i/MSM (primordial abundances of the singlet leptons) can be generated during inflation 
or reheating of the Universe. For this end we added to the SM and to the i^MSM higher 
dimensional operators, breaking these conservation laws, and suppressed by the scale QMp. 
The bounds on parameter Q were derived from requirement that the Higgs inflation is not 
spoiled. The constraints on Q allowed to estimate the B — L excess and the abundances 
of singlet fermions after reheating. We found that the contribution of higher dimensional 
operators, bounded in this way, to the production of Warm Dark Matter (WDM) and to 
baryogenesis is negligible. We concluded, therefore, that the production of sterile neutrino 
WDM 2 and baryon asymmetry of the Universe must be a low-temperature phenomenon, 
having nothing to do with inflation and reheating. 

The aim of the present paper is to revisit this analysis in connection with new results 
obtained recently in the investigations of consistency of Higgs inflation. We present below 
a short account of this progress, to motivate our present study. 

For successful inflation yielding primordial density perturbations with the amplitude 
matched to the cosmological observations [5], it is required at the tree level that the non- 
minimal coupling £ must be rather large, 



where A is the Higgs field self-coupling constant. This relation is modified by loop correc- 
tions, see [6, 7, 8, 9] for numerical results and [10] for further discussion. 3 

The requirement £ 3> 1 raised doubts in the validity of the Higgs inflation [11, 12]. The 
estimate of the high energy scattering amplitudes in the low energy electroweak vacuum 
showed that the theory enters into a strong-coupling regime (the tree unitarity is broken) 
at energies E well below the Planck scale, E ~ Mp/£. Since this energy is much smaller 
than the values of the Higgs field during inflation, H ~ Mp/\/£, the authors of the papers 
[11, 12] concluded that the Higgs inflation is an "unnatural" phenomenon, as it intrinsically 
requires fine-tuning of the Higgs potential at large values of the Higgs field. 

However, the analysis [10] (see also [13]) revealed that the theory remains in a weak 
coupling regime below certain cutoff scale A(h), which in general depends on the Higgs field 
and is in fact above the inflationary scale. This makes the Higgs inflation self-consistent and 
"natural" [10, 13]. However, the phenomenology and predictivity of the model depend on 
the properties of the UV-complete theory at high energies [10]. Here we analyze the post- 

2 On the contrary, a sterile neutrino with the mass exceeding 100 keV (a Cold Dark Matter candidate) 
can be created during the reheating stage of the Universe in the necessary amount. 

3 Everywhere in this paper, except for the Section 3.1, we will need the values of the constants at the 
inflationary scale. Then, this relation is hardly modified compared to the tree level result, and the only 
subtlety is the connection of the Higgs self-coupling at high scale with the physical Higgs mass. Under the 
assumptions, formulated in ref. [10], this relation is quite precisely given by the SM renormalization group 
running, see Fig. 5 of ref. [7]. 




(1.2) 
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inflationary phenomenology using the effective field theory parametrization of the effects 
which might come from the high energy theory. 

One can distinguish several cutoff scales corresponding to violation of perturbative 
unitarity in the different sectors of the model. In the gravity-scalar sector the cutoff is 4 



Ag-s (h) ~ < 



f Mil f 0I h< Me. 
UlI for Me. < h < Mil fn) 



The cutoff for purely gravitational interactions (graviton scattering off graviton) is just the 
effective Planck scale (see eq. (1.1)) 

Aknck^P + ^ 2 - (1-4) 

The cutoff associated with SM gauge interactions at h > is 

Agaugc — h ■ (1-5) 
and coincides with at smaller fields h. 




Figure 1: Schematic depiction of the cut-offs (1.3), (1.4), and (1.5) in the Jordan (left) and Einstein 
(right) frames. 

One of the reasons the inflation takes place while the theory remains in weakly coupling 
regime is because it does not involve large momentum transfers. The validity of the theory 
above momenta A(/i) is unclear. Thus, the question arises what happens when momentum 
transfers E are between A gauge (h) and Api anc k. Remaining within the minimal setup, 5 we 
can imagine two scenarios: 

4 We present all cutoffs in the Jordan frame, corresponding to eq. (1.1) and use the unitary gauge, 
H T = (0,h/s/2). 

° Alternatively, it is possible to attempt perturbative UV completion by introducing new states at the 
cut-off scale [14]. 
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1. The Lagrangian of gravity + SM (actually, its viable and phenomenologically com- 
plete extension such as the z^MSM) remains valid up to the Planck scale (without 
any higher order operators suppressed by lower scales), but the system at some h 
and large momenta A gauge (h) < E < Api anc k is in the strong coupling regime and 
requires non-perturbative methods for its description. Essentially this point of view 
was taken in our previous studies of inflation and reheating of the Universe [1, 2, 6, 7]. 

2. The consistency of the theory in this range of energies requires introduction (below 
the cut-off scale) of higher order operators, suppressed by a Higgs- dependent cutoff 
A(h). 

The present paper is devoted to the study of phenomenology of the second scenario 
in the SM and z/MSM, as the first one was considered in [2]. The fact that the cutoff is 
considerably smaller than the Planck scale may induce a number of interesting effects at 
low energies, such as proton decay or neutrino mass generation. In the early Universe, 
the baryon and lepton number violating reactions are greatly enhanced in comparison 
with the first scenario, which can potentially change baryogenesis or leptogenesis. The 
sterile neutrino WDM primordial production at the reheating stage after inflation may 
rise considerably in the second scenario. At the same time, even with the small, but 
Higgs- dependent cutoff, the inflation is not modified. 

The paper is organized as follows. In Section 2 we discuss higher-dimensional operators 
which could emerge in the model with Higgs-dependent cutoff (1.3) and explain why they do 
not spoil the Higgs inflation. The case of the SM is considered in Section 3. We study low- 
energy phenomenology and reconsider the processes in post-inflationary Universe. We show 
that higher-dimensional operators may produce lepton asymmetry and hence be responsible 
for the matter-over-antimatter domination in the present Universe. We explore the possible 
role of higher-dimensional operator in reheating. Section 4 is entirely devoted to the i/MSM 
and we address here all the questions we do in Section 3. We show, in addition, that the 
model new fields — sterile neutrinos — might form warm dark matter, and describe new 
options for reheating and leptogenesis. We summarize our results in Section 5, concluding 
with discussion of UV-completion of the model with regard to cosmology. 



2. Prerequisite discussion: Higgs-inflation framework and higher dimen- 
sional operators 

We start with introduction of relevant framework and short review on the relevant facts 
about the Universe evolution in the model with the non-minimally coupled Higgs inflation 
following refs. [1, 2]. 

From now on all the computations in the paper we are going to make in the Einstein 
frame, which is connected to the Jordan frame (1.1) by the conformal transformation 

Qiiv -> = ^ 2 9iiv , ^ 2 {h) = —2 • (2.1) 
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This change induces modification of the kinetic term of the scalar field, which can be 
returned to canonical form by the field redefinition 



h for h < X CI = ■v/|4 E , 

X*< nr... . . . V £j_ (2-2) 



3 fl/f„ 1„„ n2(U\ f^v h \ Y — . 1 2 M P 

cr \/ 3 £ 



M P log ^(/i) for /i>X c 



The action (in the Einstein frame) then becomes an action of the usual minimally coupled 
to gravity scalar field x with the potential 

The rest of the action (after proper rescaling of the fermion fields) looks like the usual SM 
action in the unitary gauge with proper rescaling of the mass terms by O [h (x)]- I.e. the 
usual SM mass terms for fermion ^ and vector boson get the general form of 

_ Yh( X ) _ g 2 h\ X ) A Atl , 4 , 

^Yukawa ^ W , ^gaugc - ^ ^ ^ ^ ' V A > 

with Y and g being the corresponding Yukawa and gauge coupling constants. 

The potential (2.3) provides for slow roll inflation in the region h > Mp/vf (x > 
Mp), and gives rise to primordial density perturbations of observed amplitude for model 
parameters obeying (1.2). For small field values h < Mp/^ and equally small energies the 
model behaves as the usual SM: this is the usual thermal evolution of the Universe. 

In the intermediate region Mp /£ < h < Mp / the scalar potential is approximately 
quadratic. Oscillations of the scalar field with frequency 

[X Mp , . 

"-VaX (2 ' 5) 

and amplitude in the intermediate region above reheat the Universe. In the absence of the 
higher dimensional operators suppressed by Higgs-dependent cutoff the reheating happens 
via production of gauge and Higgs bosons by the oscillating field background as described 
in [2, 3], and the reheating temperature T r is confined in the range 



A 



1/4 



1.4 x 10 _5 M P < T r < 4.5 x 1(T 5 (^^J M P , (2.6) 

where the relation (1.2) is put in. 

The overall conformal factor Q (/t) becomes relevant only for h > Mp/y 7 ^, which 
corresponds to the inflationary regime. This means, that after inflation (starting from the 
post-inflationary reheating) the value of the cut-off in the Jordan and Einstein frames are 
the same. As far as the current work deals mostly with the post inflationary processes, 
the only change we need to make to get to the Einstein frame with canonically normalized 
fields is the variable change (2.2). 

Let us proceed with the discussion on higher dimensional operators. As far as the full 
UV complete theory for the non-minimally coupled Higgs inflation is not yet developed, 
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one can not make any definite predictions. Meanwhile we can explore what happens if 
various higher dimensional operators are added to the model. 

We will analyze the effect of the following higher dimensional terms 




(2.7b) 



(2.7a) 



(2.7c) 



where L a are SM leptonic doublets, a = 1,2,3, N stands for right handed sterile neutrinos 
potentially present in the model, H a = e a bH£, a, b = 1, 2; here A depends on the background 
(Higgs) field value, A = A (h). 

In [2] the limits on these operators were obtained in the assumption that the cut-off 
scale is constant. Then these operators either directly change the inflationary potential, 
like the Higgs potential (2.7a), either generate the contributions to the potential by loop 
corrections. For example, from the first two operators in (2.7) contributions in the Jordan 
frame are h 6 /A 2 and /i 8 /A 4 log(m 2 (fr)/^ 2 ). Both terms grow much faster than h 4 and 
can not be compensated by the fourth power of the conformal factor, thus leading to 
non-flat potential and bad inflationary properties. Numerical estimate [2] shows, that the 
minimal value of constant cut-off, consistent with inflationary predictions, is A const ~ Mp. 
However, all the choices (1.3), (1.4), and (1.5) grow proportionally to h in the inflationary 
regime, and the dangerous terms in the potential behave just as h 4 for inflation. Then 
after conformal transformation the shape of the inflationary potential (2.3) remains the 
same, with modified constants. This means, that the inflationary predictions for the initial 
density perturbations remain unchanged, and the only modification is a correction to the 
normalization (1.2) of the non- minimal coupling. 

Perturbativity up to energy scale A implies that the numerical coefficients a, (3 are 
of order one or less. With the lack of knowledge about underlying theory, this is the 
only requirement on higher dimensional terms (2.7). However, one can distinguish three 
different situations. 

First, it could be that still unknown dynamics at the scale A <C Mp does not violate it- 
self any global symmetries of the SM. Then, numerical coefficients in front of corresponding 
symmetry-violating operators are suppressed, so that the effective energy scale entering, 
say, baryon number violating operators in (2.7) is the Plank mass or even higher. The 
opposite case is when all gauge invariant operators enter (2.7) with numerical coefficients 
of order one. Finally, the intermediate situation corresponds to some hierarchy in these 
coefficients. As we observe below, this third situation is the most interesting, as such 
models exhibit the richest phenomenology. In particular, we will discuss in due course the 
case where at the scale A higher order operators could dominate over similar ones of lower 
order. In particular, for the operators of the form 



5C T = y T L T HE T + (3 y L T HE T 



(2.8) 



A 2 



+ ••• 
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one can fancy the hierarchy 

1 ~ P v > Vt ~ • (2.9) 

For small H this gives nothing new. However, during reheating we can get H ~ A, so the 
second term in (2.8) can overcome the first one at large values of the Higgs field, and make 
a significant contribution compared to the lowest order estimate. Note, that hierarchy 
(2.9) does not generally imply strong coupling below the Plank scale. To simplify further 
discussion we suppose that all similar higher order operators are negligible. 

Note in passing that in this work we treat a particular set of couplings /?, 06, ■ ■ ■ from 
pure phenomenological point of view without any relation to what happens to the operators 
(2.7) above the scale A and without any discussion of matching conditions at critical values 
of the Higgs field (1.3). Having in mind a possible UV-completion like asymptotic safety 
[15] we also disregard all dangerous contributions naively expected in the theory (2.7) with 
a formal high-energy cutoff at the scale A, such as quadratic in A corrections to the Higgs 
boson mass and others. 



3. Higher dimensional operators in the SM 

3.1 Role of the non-renormalizable operators in low energy (present day) physics 

At present time the Higgs field dependence of the cutoff in the (2.7) is not observable, since 
the Higgs background has relaxed to its vacuum value h = v = 246 GeV. So the effect of 
the higher dimensional operators (2.7) can be analyzed as usual, and is most important for 
the lowest values of the cutoff, that is Mp/£, see eqs. (1.3)-(1.5) and the nearby discussion. 
As far as the operators are suppressed by this scale, the experimentally interesting effects 
can be expected from the operators violating symmetries of the renormalizable SM action. 
These are the operators in (2.7b) responsible for the violation of lepton and baryon numbers. 

The lowest order non-renormalizable lepton number violating operator (first term in 
(2.7b)) provides the Majorana mass terms for active neutrinos 

^ 4A — 2*P ^ ' 

To explain the observed pattern of the active neutrino masses the scale of this lepton 
number violating operator should not be higher than 

A < 2.8 x 10 14 GeV x p L x ( 3x 10 ' 3eV ^ ^ . (3 . 2) 

V Amf tm ) 

For the allowed range of parameter £ at the electroweak scale (see [7] , Fig. 4 for discussion 
and relation to the Higgs boson mass), £ ~ 10 3 — 10 5 , we get for the scale A 

A = Mf. „ (0.3 - 35) x 10 14 GeV . (3.3) 

The neutrino masses can be explained without taking too large parameter j3, namely with 

f3 L ~ 0.1 - 10 . (3.4) 
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The next operator, important at low energies, is the baryon number violating one, as 
it leads to the proton decay. We can use the estimate of the proton decay in the standard 
SU(5) model to place a lower limit on the scale A. Indeed, for the SU(5) case the proton 
life-time is [16] (the exact form of the relevant dimension-6 operators contributing to the 
decay is given in the reference) 

Ws+ = 1.0 x 10»yea rs x (V*.)* {J^f , (3.5) 

where My = and v$ is the vev responsible for spontaneous breaking of SU (5) gauge 

symmetry. Thus, current limit on the proton decay, r p _ 5 . 7r o e + > 1.6 x 10 33 years [17], yields 
a bound on the scale t>5, which in our case corresponds to the scale A so that 

A > VPs x 10 16 GeV x ( ) ^ . (3.6) 

\ 1.6 x lO w years / 

The discrepancy between the limits (3.2) and (3.6) is about 40 (for (3b = (3l = 1) and 
can be reconciled only with a particular hierarchy between the dimensionless coefficients: 

0L > 36 x ^/Jb . (3.7) 

For example, if (3l ~ 1 within the interval (3.4), one has to chose (3b < 10 -3 . Thus, baryon 
number violating scale cannot be A = Mp/£, it is certainly higher. 

This implies, that while the neutrino masses can be explained for a wide range of 
the Higgs masses, if the lepton number is violated at the scale (1.5) or (1.3), the proton 
decay should be prevented by some means. If allowed, its rate must be suppressed by a 
higher scale, like gravity scale (1.4). The required hierarchy between (3l and \fj3~B though 
large, still may be accidental; then it suggests that proton decay has to be observed in the 
nearest future experiments. Otherwise, the hierarchy may be a natural outcome of some 
unknown mechanism operating at the scale of UV-completion. Indeed, for instance, the two 
operator under discussion (2.7b) differ in a particular respect: one involves the Higgs fields, 
while another does not. A specific renormalization flow evolution of the Higgs-containing 
operators might be responsible for establishing the required hierarchy (3.7) at the scale 

3.2 Neutrino reheating via dimension-5 operators 

It was shown in [2], that if the suppression scale equals the Planck mass, then the particle 
production by the higher dimensional operators is negligible, and hence their impact on the 
reheating. However, this is not immediately obvious with the Higgs-field dependent cut-off. 
Let us analyze the contribution to reheating of the lepton number violating operator in 
(2.7a). 

Scalar gravity sector cutoff. We start from the operator (3.1), and omit the numerical 
coefficients F a p ~ 1 in further estimates. During the preheating, while the background 
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field evolves according to eq. (3.4) of [2], we have the following effective time-dependent 
Majorana type mass source for the active neutrinos 



&tJ\u , for x(t) > Jfw , 



&,/AI Y 2 f or y (V) < 



While the number of generated fermions of given 3-momentum k is small, <C 1, the 
effect of this can be treated by first order perturbation theory. The calculation described 
in Appendix A yields for the fermion density saturated by the moment of reheating 



Without a hierarchy in F a p the total number density of all three species of active 
neutrinos is roughly trice the amount of (3.9). To check, whether this leads to preheating, 
we should compare the energy density in neutrinos with that in the inflaton oscillating 
with amplitude X, see Appendix A for details. The equality between these two culminates 
the reheating, and the corresponding reheating temperature would be 

10 Ay/* fj L M P 



Though it is parametrically larger than the temperature of reheating due to renormalizable 
operators [2], j^ enorm oc but numerically, for (3i ~ 1, it is of the same order, cf. (2.6) 

T r ~ 1.7 x p L x 10 14 GeV ~ T T r cnorm ~ 10 14 GeV . 

Hence, higher order operators do not significantly enhance reheating. 

We check in Appendix A that the Pauli blocking is relevant for /3l ~ 1 if the reheating 
temperature is not at the lower end of the allowed range (2.6). This diminishes the possible 
contribution of higher dimensional operators to reheating: they are never really relevant 
to this process. 

Gauge sector cutoff. One could think, that with the lower cut-off (1.5) the contribu- 
tion to reheating becomes significant. However, with more active generation Pauli blocking 
effects would suppress the process much earlier, so the contribution of the generated neu- 
trinos to the energy balance during reheating can not really exceed calculations in the 
previous section. Of course, as before, one would expect that still existence of additional 
reheating channel somewhat favors the higher end of the reheating temperature interval 
(2.6). 

Finally, we conclude, that the estimate (2.6) is still a valid one, with additional pref- 
erence for the higher temperature part of the allowed region. 
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3.3 Lepton asymmetry in the Early Universe from 5-dimensional operator 

The next interesting effect which can appear from the lepton violating operators in (2.7) 
is the generation of the baryon asymmetry via leptogenesis. 

One can ask whether it is possible to generate a reasonably large lepton asymmetry 
Arii = til — n>i with the nonrenormalizable operator (2.7b) at the post-inflationary pre- 
heating stage. At a given time the amount of asymmetry is determined by averaging of 
lepton number operator Ql, which evolves in accordance with the Heisenberg equation, 

i -7: Ql = Hi n t, Ql 
at L 

Here the interaction Hamiltonian comprises both SM interactions and all nonrenor- 
malizable terms (2.7). From these equations one finds in a usual way (see, e.g. [18]), that 
at a given time the production rate of lepton asymmetry, dAriL/dt is proportional to the 
imaginary part of the minimal nontrivial trace of the SM Yukawa couplings Y a a = Y a ■ 5 a p 
entering 

C Y = —Y a L a HE a + h.c. , (3.11) 
and to the couplings F a p entering (2.7b). This trace reads 

Pi Tr (FF^FYYF^YY\ , (3.12) 

and, since the Yukawa constants Y a p can be chosen to be diagonal, the dominant term is 

/3iYfF 3 pF*pF a3 FZ 3 , (3.13) 

where I3 = y T obviously provides with the largest contribution. For simplicity we adopt the 
unitary gauge while making an estimate of the generated asymmetry. The Majorana terms 
(2.7b) connect only neutrinos, while the Dirac terms (3.11), instead, connect electrically 
charged leptons. Thus, to get the term (3.12) from a corresponding fermionic loop one 
has to insert the electroweak vertices with VF-bosons everywhere between the vertices 
proportional to Y and F (overall four electroweak vertices). These TV-bosons can be either 
virtual or external. The former case implies 3-loop contributions, while the latter case 
actually describes scattering of VF-bosons produced in the early Universe via Higgs field 
oscillations as explained in ref. [2]. 

Let us consider the 3-loop contribution, which is numerically suppressed by the loop 
factor 

During the reheating the relevant energies for the generation asymmetry are given by the 
Higgs field oscillation frequency co, which is much smaller than the gauge boson mass at 
this epoch, see (2.4), (2.5) 



/A M P M P \ X \ 



3 i a v 
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Thus, W-bosons can be actually integrated out which gives the effective Fermi lagrangian 
with the following combination playing the role of the squared Fermi constant (see eq. (3.8) 
in [2]) 

9 A e i 



- Mix 2 ' ^ ■ (3 - 15) 
The couplings (2.7b) and (3.11) yield the Higgs-dependent factor (see eq. (3.9) in [2]) 

X 2 M| f x M P \ l ../ft 2 



e v «a i " UJ ■ < 3 - 16 > 

Collecting the factors (3.14), (3.13), (3.15), and (3.16) we get for the asymmetry generated 
during one oscillation of the Higgs field 

Let us now estimate the result for the possible choices of the high-energy cutoff scale. 

Gauge cutoff. The maximal effect is expected from the lowest available cutoff, the gauge 
one (1.5). In this case the ratio h/A goes to one, and asymmetry generation rate scales 
as An^ oc h 4 oc \ 2 - The amplitude X of the oscillations of the field \ is inversely 
proportional to the cosmological time. This means, that integration of the asymmetry 
over a period of matter dominated expansion corresponds to multiplication by the time 
of expansion (or number of oscillations). Recall, we found in Section 3.2 that reheating 
develops here almost as it does without any higher dimensional operators. For reheating 
culminating at the amplitude of the Higgs field oscillations X r (see [2]) the number of 
field oscillations during reheating is N ~ £X cr /X r , and the entropy after reheating is 
s ~ (?y 4 A 3//4 (X r X cr ) 3 / 2 . Collecting this we get for the lepton asymmetry 

For the estimated in [2] range of reheating, caused by lower order operators, eq. (3.17) 
gives the asymmetry to be in the range 

4 1 0.25V'" .. ,„-U „ t ^ a „ „4 I Vr \ 4 /0.25V' 2 „ , „ ,n-15 



or, inserting relation (1.2), 

where the smaller value corresponds to the larger reheating temperature. One can see, 
that if /3l ~ 1 this may be just on the very edge of reasonable value. In particular, for 
the ranges (3.18), (3.4) one can both explain the active neutrino masses and generate a 
sufficient amount of lepton asymmetry at preheating stage. However, the result is rather 



- 11 - 



sensitive to the higher dimensional operators in the SM lepton sector (2.8). As far as during 
reheating h ~ A, all the higher order operators in (2.8) are not automatically suppressed 
compared to the leading one. This means, that just the largest one contributes (assuming 
that the whole sum still behaves nicely). Thus, for the example of hierarchy (2.9) the result 
(3.17) is enhanced by a factor (Py/yr) 4 leading to increase in the net lepton asymmetry 
up to 8 orders of magnitude. In this case quite large amount of lepton asymmetry can 
be produced. Note, that the higher order operators enhancing the lepton asymmetry 
generation are obviously irrelevant for lepton number violating processes at low energies 
discussed in Section 3.1 and hence phenomenologically viable. 

Gravity-scalar sector cutoff. Otherwise, if we adopt the gravity-scalar sector cutoff 
(1.3), which is the smallest one available, then the ratio h 2 /A goes to Mp/£. We are left 
with An^ oc const, so in the expanding Universe the asymmetry is saturated at lower 
times, that is at the reheating. Integration of the asymmetry over a period of matter 
dominated expansion corresponds 6 to multiplication of An^ by the total number of oscil- 
lations. Then, for reheating which happened at the amplitude of the Higgs field oscillations 
X r (see [2]), it took approximately N ~ £X CI /X r oscillations and lead to the entropy after 
reheating s ~ g , * //4 A 3/ ' 4 (A r A cr ) 3 / 2 . Collecting all the relevant factors we get for the lepton 
asymmetry 

For the estimated in [2] range of reheating, caused by lower order operators, this means 
asymmetry in the range 



0.25\ 15/4 n 17 _ . „ 4 /0.25\ 5/2 
I ^ n <7 w 1 n— i't^A-^fl 4 / l 



ft [—) x 2-7 x lO-^e < A L < ft x 10" i4 e , 

or adopting (1.2) 

f3l 13/4 x 6.3 x 10- 13 < A L < ft f^j 2 x 2.4 x 10' 10 (3.20) 

where smaller value corresponds to larger reheating temperature. One can see, that if 
Pl ~ 1, this is too small asymmetry. However, the result is rather sensitive to the higher 
dimensional operators in the SM lepton sector (2.8). By the end of reheating we have 
H^H/A 2 ~ 3X CI /4X r and all the operators in (2.8) contribute proportional to their cou- 
pling constants. Thus, dimension-6 operator (2.8) with hierarchy (2.9) still leads to en- 
hancement, but not that large as with the gauge cutoff, and decreasing with higher re- 
heating temperatures. The result (3.17) is enhanced by a factor (3f3 y X cr / '4y T X r ) 4 leading 
to increase in the net lepton asymmetry: lower and upper limits of (3.20) grow by factors 
of 15 and 2 x 10 5 , respectively. We should keep in mind, that the higher value of the 



With accuracy of a factor of order one, which depends also on the treatment of field-dependent cutoff 
in the calculation. 
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reheating temperature is more realistic (especially with account of additional contributions 
from neutrino reheating described in Section 3.2), and at high reheating temperatures even 
this enhancement does not lead to viable leptogenesis. So, generation of sufficient baryon 
asymmetry with scalar-gravity cutoff is possible, but may be problematic, unless (3l takes 
highest values in the interesting interval (3.4). 

Finally, at later stages of the Universe evolution the first term in (2.7b) gives rise to 
lepton violating scattering in plasma [19], which might wash out the lepton asymmetry 
generated at preheating. To clear up the situation we adopt formulas from ref. [20] for 
our case of v r playing the major role in lepton asymmetry generation. With parameters in 
front of the first term in (2.7b) yielding the correct values of masses and mixing angles in 
neutrino sector we found that the lepton violating scatterings are slower than the Universe 
expansion rate at T > 10 14 GeV and hence are out of equilibrium. This temperature is 
about the reheating temperature of the model (2.6), and the washing out processes (if any) 
are not efficient. At later times the sphaleron processes come into equilibrium and transfer 
the lepton asymmetry generated at preheating to the net baryon asymmetry [21] As ~ 
A^/3. From the estimates above, (3.18), (3.20) we conclude that the baryon asymmetry 
of the Universe can be explained with higher dimensional operators, if suppressed by the 
gauge cut-off (1.5), and with some hierarchy like (2.9) this conclusion becomes certain. 
In the case of gravity-scalar cutoff (1.3) the hierarchy like (2.9) is definitely required in 
order to produce the sufficiently large amount of the baryon asymmetry. Moreover, the 
reheating temperature in this case has not to be too high, otherwise another mechanism 
of baryogenesis has to be implemented in the model to explain the observed matter-over- 
antimatter asymmetry of the Universe. 

Note in passing that in principle additional operators in quark sector can also con- 
tribute in a similar way to the generation of the baryon asymmetry, but these operators 
are proportional to higher powers of H^H/ A 2 and the generic analysis is rather involved. 

We conclude this Section by the statement that with higher order operators in Higgs- 
inflation there are possibilities to explain both the neutrino oscillations and baryon asym- 
metry of the Universe. 

4. Higher dimensional operators in the i/MSM and dark matter sterile 
neutrinos 

DM problem needs introduction of some new particles in the SM. A minimal example of 
such physics is z^MSM [22, 4], the Standard Model extension with three sterile neutrinos 
Nj, I = 1,2, 3, lighter than electroweak bosons. To the SM Lagrangian the following terms 
are added 

5L „MSM = iN ld ^ Nl + I MjNfNj - f aI L a NjH + h.c. (4.1) 

The lightest sterile neutrino, say mostly N\, comprises dark matter while the two heavier 
sterile neutrinos due to nonzero Yukawa couplings f a j <C 1 provide active neutrino with 
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masses. Sterile neutrinos can contribute to the dark matter and lepton asymmetry gen- 
eration. One of the important difficulties in the plain z^MSM setup is the generation of 
correct abundance of the DM sterile neutrino (see [4] for a review), which requires very 
severe fine tuning in sterile neutrino sector [23] or a special source of dark matter neutrino 
in the early Universe, e.g. decaying inflaton [24, 25]. We consider in this Section whether 
the higher order operators (2.7) can change the situation. 

The sterile neutrino has renormalizable Yukawa interactions (4.1) with the fields of 
the SM. These interactions inevitably give rise to the decays of N such as N — > 3v and 
N — > jv, and to its production in the early Universe. The dark matter candidate in the 
form of sterile neutrino was introduced in [26] (for a review see [4] and references therein). 
The analysis of the non-resonant production of neutrinos and the X-ray constrains puts an 
upper limit on its mass, M < 4 keV, which happens to be in conflict with the lower limits 
on from Ly-a analysis M > 8 keV. 7 

We will see, that the higher-dimensional operators can change the situation and lead 
to efficient production of this DM candidate, consistent with other constraints. Still, even 
in this case the mass of sterile neutrino should be relatively small, as is argued below. 

4.1 Stability at cosmological time-scales 

The assumption that all higher dimensional operators are included is very powerful as it 
allows to conclude generally that any DM candidates have to be unstable (though very 
long-lived) and relatively light, so that the best strategy for their search is searches for 
emerging in dark matter decay photons of MeV or lower energies. Considering i/MSM as 
an example, the argument above goes as follows. 

Even if all couplings f a \ in (4.1) are zero, they will be generated by five-dimensional 
operators, proportional to bj Ja in (2.7). The use of equations of motion for N leads to an 
estimate of effective Yukawa coupling 

Mi 

/al~ky- (4-2) 

So, N\ is unstable due to higher dimensional operators, with the partial width of the 
order 

9 6| aG 2 F v 2 M\ 

Tn ^ ~ ~512^ AT • (43) 

Taking a constraint from EGRET data t 1v > 10 27 s[30], we obtain that Mi < 200 MeV 
for bL a ~ 1 and A = Mp. Similar bounds follow from five-dimensional operators like 
the second term in (2.7c) where derivatives act on lepton doublets. With smaller scale 
A = Mp/£ instead of Mp in this term (or larger than (4.2) Yukawa couplings), the upper 
limit on sterile neutrino mass becomes stronger, 

Mi < 4 MeV . 



7 The resonant production of DM sterile neutrinos [27, 28], which occurs in the presence of lepton 
asymmetry, escapes these constraints [29]. 
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It is also worth mentioning that light dark matter sterile neutrino can not help to explain 
the active neutrino oscillations, as it contributes only a tiny amount to sterile neutrino mass 
via seesaw mechanism in (4.1); to gain active neutrino masses, heavier unstable neutrinos 
are required, which however can be searched for in particle physics experiments [31]. 

4.2 Dark matter production 

The sterile neutrino N% can be generated from the dimension 5 operator in (2.7) during 
the hot stage and at the preheating stage. 

Production at hot stage. The coupling, proportional to (3n in (2.7), produces sterile 
neutrinos in the primordial plasma at hot stage of the Universe evolution. If stable at 
cosmological time scales this sterile neutrino is a viable candidate for dark matter. If it 
never came to equilibrium, its mass is limited from above at a given /3tv/A by imposing 
$7 at < Qdm. This limit can be read off from eq. (5.20) of ref. [2] after multiplying the l.h.s. 
of that equation by 4A 2 /Mp and choosing T r = A and f3 = (3^. This gives 

Equality in the equation above saturates if sterile neutrinos comprise right amount of the 
dark matter. 

We observe that with A ~ Mp and /?at ~ 1 the dark matter neutrino has mass in 
sub-MeV region. The lower the scale A, the lighter the dark matter neutrino. With low 
energy cutoff A = Mp/£ one obtains 

/n nQ\ 2 
M N < 10 x — keV . 



Neutrino average momentum is close to the thermal average momentum in the primordial 
plasma, so the constraints from Ly-a discussed in [4] are relevant in our case either. Thus, 
for the cutoff A = Mp/£ the viable dark matter implies /3jv < 0.03, sterile neutrino might 
form warm dark matter. If f3j^ ~ 1 and A = Mp/£, stable at cosmological time scales 
sterile neutrinos contribute to the hot dark matter component, and hence can not solve the 
dark matter problem. With /3jy ~ 1, A = Mp the sterile neutrinos are cold and can solve 
the problem. Note, that a particular value of dark matter mass Mn should be corrected 
if higher order operators change somewhat the reheating in the model, as we discuss in 
Section 3.2. 

However, we will see now, that at the preheating stage the neutrino production is much 
more active. 

Dark Matter generation during preheating with dimension-5 operator sup- 
pressed by A(/i). Let us consider the sterile neutrino production by the same dimension 
5 operator as that in previous Section, suppressed by the Higgs-dependent scale, corre- 
sponding to the gravity-scalar cutoff (1.3), 

^=^^ N =Twr N - (iA) 
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Calculation proceeds along the same lines as in Section 3.2. The sterile neutrino density 
at late times is given by eq. (A. 7) with obvious replacement fit — > (3/y. 

Dividing this density by the entropy at reheating we get for sterile neutrino abundance 

An{X) ~ ^M^U^^Tf 



90 q2 V\ 2 3 / 2 M%_. - in _ 4fl2 H^lMl 



(seemingly large £ suppression is not there, as it is compensated by two extra powers of 
Mp/T r ). Thus, we get for Q, DM = 0.223 

n N (3% M N / A 1 Mp 2 M N /025 fl.l x lO^GeV^ 3 



n DM 2.4 keV V 0.25 £ 2 T r 3 ^ N 1.3 x 10~ 4 keV V A \ T r ) ' (46) 

We can see, that if the constant /3yv > 10 -2 , the sterile neutrino DM should be lighter than 
keV, which violates the structure formation bounds. Thus, if f3j^ ~ 1 it is a problem to 
generate DM sterile neutrino with the higher dimensional operators; however with hierarchy 
Pn 10 -2 it is still possible. Obviously, the case of gauge cutoff leads to even larger 
production. 

The case of the Planck mass cutoff (1.4) coincides with the analysis in [2], which also 
allows successful generation of the DM density 



n N 2 M N /0.25A /l.l x 10 14 GeV\ 
JW ~ 7V 10 4 keV {^) { T r ) ' 



(4.7) 



In this case the sterile neutrino should be rather heavy, and additional care is needed to 
assure its stability (see previous section). 



5. Summary and conclusions 

We have analyzed various effects that can appear from the higher order operators in Higgs 
inflationary theory. In the previous analysis [2] we used additional operators suppressed 
by a constant (field independent) cutoff, which lead to the conclusion that either the 
inflationary properties of the theory are spoiled, or no significant generation of any chemical 
potentials is possible at preheating due to these operators. 

In this paper we analyzed a wider class of operators, suppressed by the Higgs dependent 
cut-off, which is essential to sustain good properties of the theory [10]. The key observation 
here is that now the cut-off can be lowered significantly at present time and at reheating 
without spoiling the inflationary properties. As a result, with addition of properly chosen 
higher dimensional operators one can explain many effects. With proper tuning of the 
coefficients in front of these operators, one can explain neutrino masses, obtain proton 
decay, generate baryon asymmetry and sterile neutrino DM (within i/MSM framework). 
However, all this is only possible with rather specific choices of the hierarchies between the 
higher dimensional operators, and a generic prediction seems to be impossible. 
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A. Generation of fermions by varying mass term, scalar gravity sector 
cutoff case 



Let us analyze the generation of fermions by the operators (2.7b) or (2.7c). In the presence 
of the time varying Higgs field background this leads to the Hamiltonian for the fermions 

flint = J d 3 xm(t) m 

with the time varying mass (3.8). The analysis of the generation of the fermions, while 
the amount of generated particles is small (occupation numbers much smaller than unity) 
can be done perturbatively and closely follows appendix F of [2]). We treat the constant 
mass mo = f3£y~^U)/4 as a part of the free Hamiltonian. Then the number density at the 
moment t 



n k (t) 



(m(t') - m )e 2ikot ' dt' I (m(t") - m )e 



-2ik t" ^£1 



and total particle number 



n(t) 



d 3 k 



n k (t) . 



Integration is the same as in [2], leading to 



nk(ti) = 

where the constant is slightly changed: 

1 



1 



16 g j n 2 ( Trfco 



(A.l) 



L 



^ L V 3 Auj 



sin 



( — + 1 I 



1 



ko — to 



sin 



4tt 



A 1 



(A.2) 



where we removed 2k$ from the denominator and changed 3XX — > Pl\J 34^^ (compare 
(3.8) and (C.l) of [2]); the ratio uj/Uq should always stay, because it is the ratio of the 
source frequency and the energy ko in the oscillating exponent. The parameter e is defined 
as 



sin( 27 r e W||. 



(A.3) 



At large times the ratio leads to the sum of delta functions over momenta, and integration 
over momenta gives 

2 



t 



Pi 



48" 



xA (2^53 E v 1 - h " (m sin [An {l + 1} e] - rh sin [47r {l ~ 1} e] 



(A.4) 
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Summation for small e can be performed using the integral (cf. appendix E of [2]) 

f°° *du o 27T , 

S 3 = (47re) 3 — r(2iicosu-2sinn) = — (47re) 3 . (A.5) 
Jo u 3 

This finally yields the generation rate (we reintroduced also the expansion of the Universe 
here) 

Integrating this during the matter dominated expansion (3.4) from [2] we obtain the neu- 
trino density (saturated at late times) 



"-^A^- (A ' 7) 

Note, that this does not depend on time, because rate (A. 3) is proportional to X oc t~ l . 
This formula should be compared to njv(X) = ■£^XlX from Section 5.2.2 of [2]. 

Thus, the total active neutrino density is trice the amount of (A. 7). To check, whether 
this leads to preheating, we should compare the energy density in the neutrinos to those 
in the inflaton. To do this, we need energy density instead of the number density. Just 
doing the computation using occupation numbers (2.3) leads to divergent result. This is 
caused by discontinuous approximation to mass formula (3.8). However, for the estimate 
purposes we can just consider, that the average energy of created neutrinos is still given 
by the value in appendix C of [2] , i.e. 8 

E~ — = ±\fcx. (A.8) 
4vre 2 V 3 y ' 



Then we have for the energy density 

p ~ ^&%X . (A.9) 
Equating this with the inflaton energy density (3.7) in [2] we get 



^X c 3 r X = \XIX> . (A.10) 



8 One can avoid the divergence in the calculation of energy by smoothing the effective mass (3.8), say, 
as follows 



_§_L_ 3_ 
,,s 2^2 V 2A 

m - m(t) = -. — — — r- 

cosh (^JsX ■ t) 



Then, instead of (2.3) we get 



_ QjrPhl sin 2 (21^ 



32A2X2 cosh 2 (^^o) sin2 (^) 
Here we obviously get an exponential energy cut-off with typical scale (A.8). 
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Thus, the reheating (half energy in the relativistic active neutrinos) happens at 



X ~ 2lt^ cr ■ (A - U) 
This is parametrically earlier, than (3.10) in [2]. Thus, reheating temperature is 



Tr-f^Y'h!*. (A.12) 
\27gJ 2ttV^ VI 

This is parametrically larger than the temperature of reheating due to renormalizable 
operators [2], 

rpienorm Mp 

but for Pl ~ 1 it gives the numerical estimate of the same order, 

T r ~ 1.7 x fa x 10 14 GeV ~ T* eDOim ~ 10 14 GeV . 

Above we neglected the Pauli blocking effect for the fermions. Indeed, with estimates 
(A. 7) and (A. 8) we get for the occupation numbers at reheating 

- ID 3 x P L 

This reaches one (breakdown of perturbative approximation because of Pauli blocking 
effects) for (3l ~ 1 at the higher reheating temperature in [2] (see formula (3.21) there). 
So, the Pauli blocking effects were not really relevant for this calculation. 
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